Dynamics of the chiral phase transition at finite chemical potential 



O 

o 
o 

(N 

< 



(N 
> 

^t- 

(N 
l> 
O 
O 
o , 

a: 
i 

Q . 
^ ■ 



X 



Wenjin Mao 

Department of Physics, Boston College, Chestnut Hill, MA 02167 

Fred Cooper and Anupam Singh 

Theoretical Division, MS B285, Los Alamos National Laboratory, Los Alamos, NM 87545 

Alan Chodos 

Department of Physics, Yale University, New Haven, CT 06520-8120. 
Present Address: American Physical Society, One Physice Ellipse, College Park, MD 20740. 

(February 1, 2008) 

We study the dynamics of the chiral phase transition at finite chemical potential in the 
Gross-Neveu model in the leading order in large-iV approximation. We consider evolutions 
starting in local thermal and chemical equilibrium in the massless unbroken phase for con- 
ditions pertaining to traversing a first or second order phase transition. We assume boost 
invariant kinematics and determine the evolution of the order parameter er, the energy den- 
sity and pressure as well as the effective temperature, chemical potential and interpolating 
number densities as a function of r. 
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The phase structure of QCD at non-zero temperature 
and baryon density is important for the physics of neu- 
tron stars and relativistic heavy ion collisions. The phase 
structure for two massless quarks [1] reveals a rich struc- 
ture. At low temperature and chemical potential, the 
ground state has broken chiral symmetry. At higher 
chemical potential one finds a superconducting phase. 
The transition out of the chirally broken phase as one in- 
creases the temperature is second order at low chemical 
potential and then changes to first order as we increase 
the chemical potential [2]. 

Recently we found a simple model which has a simi- 
lar phase structure [3] to that described above, i.e. both 
chiral and superconducting transitions as well as asymp- 
totic freedom. Here we consider a special limit without a 
superconducting phase, where the model reduces to the 
Gross-Neveu (GN) model [4] whose Lagrangian is 
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where S~ 1 (x, y)[cr] = {"f^d^ + a) S(x — y). 

The phase structure of the GN model at finite temper- 
ature and chemical potential in this approximation has 
been known for a long time [3] [5] and is summarized in 
Fig. 1. 
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which is invariant under the discrete chiral group: ^ — > 
75 ^i. In leading order in large N the effective action is 
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FIG. 1. Phase structure at finite temperature and chemical 
potential \x. The phase below the line has (^ij)) ^ 0. 

The phase structure is determined from the renormal- 
ized effective potential 
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Here m/ is the physical mass of the fermion in the vac- 
uum sector. The tricritical point occurs at — = .608, 
— = .318. We have chosen to renormalizc the effective 
potential so its value at T = in the false vacuum a = 
is zero. In the true vacuum u = mj the energy density 
has the value e/m'j — — j^. 

Following a heavy ion collision, the ensuing plasma ex- 
pands and cools traversing the chiral phase transition. 
In hydrodynamic simulations of these collisions, a rea- 
sonable approximation is to treat the expansion as a 
1+1 dimensional boost invariant expansion [6,8] along 
the beam (z) axis. In this approximation, the fluid veloc- 
ity scales as z/t. In terms of the variables fluid rapidity 
f] = \ In (fzf) and fluid proper time t = (t 2 - z 2 ) 1/2 , 
physical quantities such as a, e become independent of 77, 
as discussed in refs. [6,8] and applied to the problem of 
disoriented chiral condensates in ref. [9]. We note that 
related noncquilibrium techniques have also been devel- 
oped in ref. [11] and applied to the problem of disoriented 
chiral condensates in ref. [12]. Although the effective 
mass a is a function solely of t, two-point correlation 
functions depend on r\ as well. 

We shall use the metric convention (—,+). In our ap- 
proximation, the dynamics are described by the Dirac 
equation with self-consistently determined mass term. 
Rescaling the fermion field, il>(x) = -^$(x), and intro- 
ducing conformal time u via r = — , we obtain 



[j°d u + 1 3 d v + a(u)] $(a;) =0, 
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where a(u) = gt = ^e u and riif = rrifT. 

Further letting g 2 — X/2N we have the gap equation 
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where we have assumed N identical = i/j. 

These equations are to be solved subject to initial con- 
ditions at r = To. It is sufficient to describe the initial 
state of the charged fermion field by the initial particle 
and anti-particle number densities, which we take to be 
Fermi-Dirac distributions described by /j,q and T . 



Expanding the fermion fields <& in terms of Fourier 
modes at fixed conformal time u, 
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the <j) k then obey 
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The superscript ± refers to positive- or negative-energy 
solutions. Introducing mode functions (f> k {u) via 



4>k( u ) 



/±(r) X ± , (0.8) 



where the momentum independent spinors are chosen 

to be the orthornomal ±1 eigenstates of we obtain 
the second order equations: 

-4-a-u>i±i < ^)f k t {u)=0, (0.9) 



du 2 



du J 



k 2 



a (u). We parameterize the positive- 
energy solutions f k in a similar manner to Eq. (3.1) of 
Ref. [10]: 
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Using eqs.(0.6, 0.8)and the definitions: (b^(k)b(q)) = 
2TrS(k-q)N + (q) and (tf(k)d(q)) = 2n5{k - q)N_(q), we 
obtain for the gap equation 

/dk 
^(l-N + (k)-N_(k))R k (u) (0.10) 

where R k (u) = 1 - 2k 2 |/+(u)| 2 . and 
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This equation is solved simultaneously with eq. (0.9). 

We take our initial state to be in local equilibrium so 
that N±(k,ti,T) = [ e ( w *( )=Ff)/T + 1 ]-i where UJk (o)=E 
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Since we start our simulation in 



= ^k 2 + a 2 {0) 
the unbroken mode, ct(0) = 0. We choose the initial r = 
^- and measure the proper time in these units. We use 
adiabatic initial conditions on the mode functions /, i.e. 



fk(0) 



/+(0) = -^ fc / fe + (0) and N 2 = [uj k (0) + 
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tr(O)] -1 . To obtain non-trivial dynamics in this mean field 
approximation at high temperatures, it is necessary to 
explicitly break the chiral symmetry by giving a a small 
initial value which we choose to be ct(0) = 10~ 3 . 

We have studied three separate starting points on the 
phase diagram of Fig. 1 in our numerical simulations. 
We determined the energy density and the pressure from 
the expectation value of the energy momentum tensor as 
described in [8]. In the 77, r coordinate system T^ v is di- 
agonal which allows us to read off the comoving pressure 
and energy density. After renormalization we obtain 
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The integrations involve a moving cutoff A = At when 
the mode functions arc truncated at physical k z = A. 
In the massless phase, one finds that the exact equation 
of state is p — e. To compare our field theory calcula- 
tion with a local equilibrium hydrodynamical model we 
assume 
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The conservation law of energy and momentum T Q/3 .p = 
0, combined with scaling law v = z/t and p = e yields 
[6] j- = (^) 2 ,^ = (^). From Eq. (0.11) we can also 
determine p(fj,,T) and e(fi,T). Assuming T/Tq = tq/t 
and /j/mo = t o/t we find that the local equilibrium ex- 
pressions for e and p evolve identically to the numeri- 
cally determined field theory evolution before the phase 
transition. (We note that in thermodynamic equilibrium 
dT/T — d(i/n [7] and so close to equilibrium we expect 
the temperature and chemical potential to have a similar 
falloff with time. In fact, a different falloff for the two 
quantities as a function of time can be viewed as a depar- 
ture from local thermal and chemical equilibrium.) With 



the same assumptions we find the distributions for N± 
plotted against k v arc independent of r. This also agrees 
with the exact evolution before the phase transition. 

We want to understand how the particle number dis- 
tributions evolve in time. In relativistic quantum me- 
chanics, particle number is not conserved. However in a 
mean field approximation one can define an interpolat- 
ing number operator which at late times becomes the 
outstate number operator. By fitting the interpolat- 
ing number densities for both fermions and antifermions 
to Fermi-Dirac distributions [13] we extract the best 
value of fi and T for that value of the proper time. 
To define the interpolating number operator we use a 
set of orthonormal mode functions y k [10] which are 
the adiabatic approximation to the exact mode func- 
■+ - "*/ Q " du : y- = v k e l I Qkdu with 
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ation and annihilation operators then become time de- 
pendent and the expansion of the quantum field becomes 

/dk 
[a(k, u)y+ («) + c t (k, u)y- (u)]e ik ^. 

This is an alternative expansion to that found in Eq.(0.6). 
The two sets of creation and annihilation operators 
are related by a Bogoliubov transformation a(k,u) = 
a k (u)b(k) +f3* k d\k); c^{k,u) =- (3 k (u)b(k)+ a%rf(k). To 
ensure that at u = the two number operators match, 
one chooses adiabatic initial conditions: <fi k = y k , so that 
Q!fe(0) = l;/3fc(0) = 0. The interpolating number operators 
for fermions and anti-fermions are defined by N + (k, u) — 
(a^(k,u)a(k,u)); N~(k,u) = (J(k,u)c(k,u)). With 
A fe = %±£ we have explicitly 
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iV±(M) = N±(k) + [1 - N+(k) - N-(k)]\p k (u)\ 2 . 

We have solved the simultaneous equations Eqs. (0.9, 
0.10) numerically. Comparing N^(k,u) with an equilib- 
rium parameterization we have determined T(k, u) and 
fi(k, u) as a function of k. When these quantities are 
independent of k v = kr this defines a time evolving tem- 
perature and chemical potential. We found that T and 
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/j, are independent of k except at high momentum before 
the chiral phase transition. 

From Fig. 2 we see that for both the 1st and 2nd order 
transitions, <t(t) shows a sharp transition during evolu- 
tion from the unbroken mode to the broken symmetry 
mode. Before the phase transition the temperature falls 
consistent with the equation of state p = e. For the 2nd 
order transition, the chemical potential follows the tem- 
perature and falls as t _1 . After the phase transition, 
there is now a mass scale m/ which leads to oscillations 
of a. For the 1st order transition the chemical poten- 
tial falls faster than r _1 . If one traverses the tricritical 
regime one finds results for [i intermediate between the 
two cases displayed. 
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FIG. 2. Evolution of T,fj, and a as a function of u. Top 
figure is for 1st order transition. Bottom figure is for 2nd 
order phase transition 

The order of the transition has a more noticeable ef- 
fect on the spectrum of particles and antiparticles. If the 
system evolves in local thermal equilibrium with a = 0, 
then when N^(k, u) is plotted vs. k n — kr it is indepen- 
dent of u. A deviation from this result is an indication of 
the system going out of equilibrium. We find because of 
the "latent heat" released during a first order transition 
that the distortion of the spectrum is greatest in that 



case, (see Fig. 3). If one traverses the tricritical regime 
one finds results intermediate between the two cases dis- 
played. 

In local equilibrium with a = 0, e = p oc t~ 2 . Simu- 
lations, shown in Fig. 4 agree with this before the phase 
transition occurs. After the phase transition we find that 
the energy density oscillates around the true broken sym- 
metry value discussed earlier, namely e = — 1/Att. These 
oscillations would be damped if we included hard scat- 
tering effects [14]. The details of this calculation as well 
as a discussion of correlation functions and the effects of 
a bare mass will be presented elsewhere. 
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FIG. 3. Evolution of N± as a function of u. Initial condi- 
tions are same as Fig. 2. The momentum displayed is k v = fcr 
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FIG. 4. Evolution of the pressure and energy density as a 
function of u .Initial conditions are same as Fig. 2. 
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